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Abstract 

Spectral decompositions for the evolution operator on an en- 
ergy shell in phase space are constructed for the free motion on 
compact 2D surfaces of constant negative curvature. Applications 
to quantum chaos and in particular to the recently proposed bal- 
listic cr-model are briefly discussed. 

1 Introduction 

The dynamics of a Hamiltonian system can either be described in terms of 
the trajectories xo —> x(t) = UtXo in the phase space, or by specifying the 
laws of evolution of a function <p(x) on the phase space: ip — > Ut<p- The 
evolution operator Ut which advances a function along the trajectories is 
defined by 

U t tp(x) = ip(U t x) = (f(x(t)). 

Quantum mechanics has a natural relation to the trajectory based 
approach through Feynman's Path Integral. In the semiclassical limit 
the path integral can be approximated by the saddle point contributions 
(which turn out to be the classical orbits) and leads to the Gutzwiller 
trace formula 0] for the Green's function of the quantum mechanical 
Hamiltonian. This is a useful tool for studying such problems as the 
quantum energy level correlations but it requires knowledge of the long 
periodic orbits in order to obtain the correlations at small energy differ- 
ences ^ The situation is particularly bad in chaotic systems, where 
the periodic orbits proliferate exponentially with length. In practice some 
uncontrolled approximations about correlations of actions for different 
periodic orbits are made to get the analytical results. 

In contrast, the flow based approach seems suitable for studying the 
behavior of chaotic systems at long times when due to the decay of cor- 
relations the dynamics becomes trivial. To use this approach in quantum 
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chaos two problems however need to be overcome. First, it is unclear how 
to relate the quantum mechanics to the evolution operator since there is 
no analogue of the Feynman path integral. Second, one needs to be able 
to calculate various properties of the evolution operator which are natu- 
rally formulated in terms of the spectral decomposition of Ut in decaying 
eigenmodes. 

The first problem has recently been addressed in j4|, || where it has 
been conjectured that the (suitably averaged) correlation functions of the 
quantum energy levels and/or quantum eigenfunctions can be generated 
from an effective action of nonlinear cr-model type involving the Liouvillc 
operator L = 4rU t (the so called "ballistic" cr-model). The inspiration 
for this approach comes from the well developed theory of weakly disor- 
dered metals, where the disorder averaged properties of an ensemble of 
macroscopically identical systems are calculated using a similar cr-model 
with a diffusion operator instead of the Liouvillian one. (See || f?J for a 
detailed discussion of the diffusive case). 

Our paper is devoted to the second problem. In particular we compute 
a generalized spectral decomposition of the evolution operator 

U t = J2e Xt \X)(X\. (1) 

A 

for a class of "model" chaotic systems namely the free motion on two 
dimensional compact surfaces with constant negative curvature. 

For mixing systems the only square integrable eigenfunctions of U t 
are the constant functions (see Appendix |A|). The eigenfunctions |A) 
entering the spectral decomposition belong to a larger space C°°(M)* 
of distributions and can be obtained from the residues of the analytical 
continuation of matrix elements of the resolvent of L (see section |^). 
The procedure was suggested by Ruelle || and was successfully used 
to study the dynamics of some chaotic maps (sec ||, |nj |ll[ fL2|). The 
eigenvalues A entering the spectral decomposition are sometimes called 
"Ruelle resonances". For our model dynamical system the resonances 
can be found using the Selberg trace formula to relate the classical and 
quantum zeta-functions This approach however does not provide 

the eigenmodes which we compute using the representation theory of 
50(2,1). 

The rest of the paper is organized as follows. We start by outlining 
the Ruelle procedure for flows concentrating on criteria for convergence 
of the resulting decomposition (section |J). We then introduce geodesic 
flows on constant negative curvature surfaces (section |^) ; summarize the 
important facts from the representation theory of 50(2, 1) and proceed 
to obtaining the spectral decomposition (section |], equation (|36|)). This 
decomposition is used to refine approximations of the decay of correlations 
for geodesic flows (section |5.l|) and relate the evolution of particle density 



on the configuration space to the Laplacian operator (section 5.2). To 



conclude we discuss the regularization of the Liouvillian operator entering 
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the ballistic c-model (section 3.2). Some technical details are relegated 
to the appendices. 



2 Decompositions of evolution operators for 
general mixing systems 

2.1 Spectral decompositions 

Consider a Hamiltonian system on a phase space T*N with coordinates 
x = (q,p) where N is a smooth compact manifold parameterised by q and 
p is the momentum. When the Hamiltonian H does not depend on time 
energy is conserved and the trajectories lie on surfaces M of constant 
energy. We study the restriction of the flow Ut to a constant energy shell 
M for some value of the energy. If dH is non-zero on M the Liouvillian 
measure dpdq induces a measure dfi on M according to dfidH = dpdq 
pHy . By Liouville's theorem this measure is preserved by the flow U t . 

The evolution operator Ut advances a function tp(x) : M — ► C along 
the flow 

U t ip(x) = ip(U t x) = e Lt ip{x) 
where for a Hamiltonian system the Liouville operator L is given by 



f d 



U t cp = {H,cp} (2) 

t=o 



with { , } being the Poisson bracket. 

The operator Ut is unitary with respect to the scalar product 



{£(x)\<p(x)) = / t(x)<p{x)dn(x) (3) 
JxeM 

since /i is preserved by Ut- A natural space for Ut is the Hilbert space 
2>2(M) of square integrable functions on M, while L preserves a smaller 
space of compactly supported infinitely differentiable functions C°°(M). 

We would like to find a spectral decomposition (Q) for the operator Ut 
in terms of its eigenvalues e xt and projectors (A| onto its eigenfunctions 
| A). This would enable us to establish the evolution of a function f(x): 



l> t |^} = ]Te At |A}(AM 



Since the operator Ut is unitary on Li2(M) it can only have eigenvalues 
which lie on the unit circle. On the other hand, in a mixing chaotic system 
(definitions and properties of mixing systems are given in appendix |Aj) all 
deviations from a constant value decay implying the existence of modes 
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corresponding to eigenvalues with modulus less than 1. The resolution of 
this apparent paradox lies in observing that the L2 (M) eigenfunctions of 
Ut do not necessarily form a basis in Li[M~) or even in C°°(M). In fact, 
for a mixing system the only square integrable eigenfunctions of Ut are the 
constant functions which have the eigenvalue 1 (appendix |A|) . In order to 
find a spectral decomposition for U t we need to extend the Hilbert space 
L^{M) to a rigged Hilbert space as described below. A similar approach 
was successfully employed for various maps (see e.g. |9[ |To|). 



2.2 Rigged Hilbert spaces 

A Hilbert space HI may be extended to the set of linear functionals on a 
suitable dense subspace §. The resulting space S* is called the rigging of 
EI over S. We denote by f[tp] the value of the linear functional / 6 S* on 
the vector tp G S. 

A vector g G H may be naturally embedded in §* as g[ip] = (tp\g), so 
we get a sequence of spaces S C EI C §* . We make the above embedding 
explicit H and denote the functional / G §>* by |/); its value on a vector 
<p G S being f[ip] = (<p\f) . We also introduce the notation (f\<p) for the 
antilinear functional f[ip]. 

We search for eigenfunctionals of the evolution operator Ut in the rig- 
ging C°°(M)* of L%(M) over C°°(M). The evolution operator is extended 
to C°°(M)* by 

(<p\tj t f) = (U-Mf) where |/) G C°°(M)* and <p G C°°(M). 

Since Ut does not preserve any scalar product in C°°(M)* its eigenvalues 
need not lie on the unit circle. 

We shall construct decompositions for the correlation function 

(Z\U t \<p) - / dtf(x)U t <p{x) 

JM 

in a subset of eigenfunctionals {|/a)} from C°°(M)*: 



{m\v)=Y, eXt ^\h){f->\v)Mfn) = e Vt \fr 1 )- (4) 

A 



In general ^ only has asymptotic meaning (see section ^3 for details), 



and will only converge when t > and £ and <p belong to a subspace T of 
C°°(M). For the free motion on compact surfaces of constant negative 
curvature we find such a subspace which is dense in the set of infinitely dif- 
ferentiable functions. Note, that the eigenfunction |/_t) with the eigen- 
value e~ At appears in (^) due to the unitarity condition (Ut) + = U-t 
which leads to another expansion for the correlation function: 



(d\U t \v) = (v\U-t\Z) =J2e- xt (£\f_j)(h\<p) (5) 
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converging when t < and (,j3£T 

Decompositions (^,||) appear naturally in connection with the corre- 
lation function (£\(L — z) -1 ^) of the resolvent of L . 

2.3 Resolvent method for calculating the evolution 
operator decompositions 

Choosing an integral representation for the resolvent of the Liouville op- 
erator converging when Re z > 

/>oc 

R-(z) = - e- zT U T dT=(L-z)- 1 (6) 
Jo 

and assuming decomposition (Q) we obtain for the correlation function of 
the resolvent 

A A Z 

where £, Lp e T. 

Conversely decomposition ([I]) , can be constructed by analytically con- 
tinuing F to the left half of the z-plane and analyzing its singularities 
and residues. The position of the poles dictate the rate of decay of the 
correlation function of the evolution operator ||. The values of A in 
decomposition (^) are called Ruelle resonances. 

Each term in decomposition is well defined for £, tp S C°°(M) 
when 

the position of the poles of F^ tV (z) 

does not depend on the choice of £ or ip. (8) 

A famous conjecture due to Ruelle || states that the poles of the resolvent 
for a mixing system do indeed satisfy this condition. 

The sum in (Q) converges if in addition the residues Res(A, F^ >v ) grow 
slowly enough when |A| — ► oo: 

lim V e At Res(A, Ft v ) = (9) 

i?,<|A|<oc 

The set T will comprise of functions in C°°(M) which satisfy condition 
@. 

To prove the convergence we consider the integral of F^, Vtt {z) — 
e zt (£\R-(z)\4>) for £ and </> in T around the contour \z\ = R for a given 
large R. In the limit R — > oo this integral converges due to (fk to 



lim / Ft v t (z)dz = 2ni VRes(A,F 4 v t ) = 2m V e At Res(A, v ) 

(10) 
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where the sum is over all the poles A other than oo. 

Since all the poles are in the left half of the z-plane the contour can 
be deformed to go along the line Re z = a for a fixed a > leading to 

F t (z)dz = e zt e- zT (£\U T \(p) dT dz 

Roz=a JRez=a J T=0 



p oo y* oo 

/ / e( a+i y^- T Uy(a\U T \y)dT 

Jt=0 Jv=-oa 



T=0 Jy=-oo 
oo 

6(t-T){t\U T \<p)dT 

T=0 

= (S\ft t \<p) fort>0 (11) 
From equation ( |l0| ) we get the decomposition into residues 

{£\U t \tp) = 27ri^e At Res(A,F) (12) 

A 

which converges absolutely for t > 0. 

The individual terms in exist for arbitrary C°°(M) but 

the series converges only for £, </? S T. The residue Res(A, F) is a linear 
functional of £ and an antilinear functional of tp. We define the operator 
K x : C°°{M) -► C°°(M)* by 

2 7 rzRes(A, J F) = <e|^ A |^) ^6C°°(^ 

so that 

A 

and -/^aIv) i s an eigenfunctional of C/j 

ft-K^) =e xt k x \ V ). 

Let {1/^)} be a basis for the eigenspace corresponding to the eigenvalue 
e At i.e. the image of C°°(M) under K\. 

K x \<p) = $>^>I/a) with ^I/a> = eA 'l/A>- 

k 

Substituting this expression into (^) we obtain 

{t\u t \<p) = YJ& k >^ 

A 

= EE^I/a)(/-aI^) ^eT,t>0 (13) 

A fc 



G 



where by virtue of the unitarity of the evolution operator the coefficients 



{ck\<fi} are given by eigenfunctionals |/ fe y) of U t with eigenvalues e xt . 



(c k \ V ) = (<p\f*J) 

Equation (|l^) reduces to (|J) in the non-degenerate case when the image 
of K\ is one dimensional. 

For arbitrary £ and <j> m C°°(M) the convergence of decomposition 

is asymptotic: 

Wt\<p)- E eAt </-»(£l/A>| <C(a)e- at ,a>0. (14) 

Re A>-a 

This inequality holds since the integral 



(f|l/t|v>- E e Xt (f-T\v)(ti\fx))e- zt dt 

Rc A > — a 

is analytic in the region Re z > —a of the z-plane. 

A similar procedure starting with the representation R+(z) for the 
resolvent converging at Re z < 

r° 

R+{z)= I e zt U t dt = {L~ z)- 1 (15) 

J — oc 

gives the decomposition (||) which converges absolutely for t < when 
the conditions (^,^) are met. 

2.4 Evolution operator at long times 

We can use the decompositions (^) and (S) to study the approach to 
equilibrium at long times for mixing dynamical systems. 

We show below that for £ and ip in T the function C(a) in ( [li] ) is 
independent of a. The behavior at long future times may be approximated 
by retaining only the terms in (Q) where A has a small negative real part. 



= E eAt </-ii^i/A)+ E eAt </-xM<£i/A> 

Rc A > — a Rc A < — a 

* E e At (/_»(£|/A> 

Rc A>-a 

= mi)(iio+ E eAt (/-A^>(ei/A>, 



0>Re A>-a 

(16) 
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where we have separated out the the contibution from the constant L2 (M) 
eigenfunction with A = 0. For a mixing system the latter is the only non- 
decaying eigenfunction. By the absolute convergence of the sum (^) for 
t > the discarded terms are bounded by 

I E e"<f-Jl\<P)(t\h)\ < e~ at £ \{f-j\vmfx)\<e- at C 

ReA<-a ReA<— a 

for some constant C and hence decay faster than e~ at . 

For long past times we work analogously from (|J) rather than (||) to 
obtain: 

(ei&b) « ^ii)(iio + E_ e-*(h\<p)(t\f-j) for * -> 

0<Rc(-A)<6 

Here the terms which decay faster than e bt as t — > —00 have been dis- 
carded. 



3 Surfaces of constant negative curvature 

We will use the resolvent method described in the previous section to 
obtain a spectral decomposition for the evolution operator of a 'model' 
chaotic system. The simplest systems which are strongly chaotic (and in 
particular mixing) are the free motion on compact 2D surfaces of constant 
negative curvature. 

The hyperbolic plane N plays the role of a universal cover for these 
surfaces. It can be embedded in Minkowski space where the metric is 

ds 2 = —dx\ + dx\ + dx\ 

as the surface satisfying the equation 

X 2 I I •^'3 — ■ 

We shall consider the free motion on compactifications of N formed by 
quoticnting it under the action of some discrete group as described below. 
We start by analyzing the free motion of a particle on the hyperbolic 

plane itself. This dynamical system has the phase space T*N and is 

2 

described by the Hamiltonian H = which is just the kinetic energy 
of the particle. The trajectories for this dynamical system are given by 
the geodesies of the surface N. For a particle with unit mass m = 1, the 
constant energy surface M with energy E = i consists of the points with 
momenta of unit modulus (and hence also unit speed). Thus the energy 
shell M is the unit cotangent bundle of the hyperbolic plane. Changing 
the energy amounts to a rescaling which leaves the trajectories unchanged 
and only alters the rate at which they are traversed. 
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The group of isometries of the hyperbolic plane G = 50(2, 1) acts 
simply transitively on the points of the energy shell. Let us take as a 
base point in M the point O = (q,p) where the position q is given by 
(0,0, 1) and the momentum p is (1,0,0). A point x on the energy shell 
may be identified with the unique element of G which takes O to x = gO. 
In this way the constant energy shell can be identified with the elements 
of G. The above construction gives a diffeomorphism from the topological 
group G to the energy shell M. 

We choose a basis for the Lie algebra g of G consisting of the matrices 





K=\0 1 (17) 



An element g £ G may be written as a product J14[ 

g(</>,T,1>)=e K +e PT e K *. (18) 

The parameters r, ef> and ip are known as the Euler angles. Note that r 
and <j) give the position in polar coordinates: 

(xi , X2 , X3) = (cosh r, sinh r cos <f), sinh r sin (f>) 

while the angle ip gives a consistent way of parameterizing the direction 
of the momentum. 

After a time t a free particle at the point O on the constant energy 
surface advances along the geodesic to which it belongs to the point h t O 
where h t = e pt . The point gO moves to ghtO; hence on identifying M 
with G the evolution corresponds to right multiplication by the element 
h t . This 1-parameter group of transformations 

9 !-> U t {g) = ght 

is the geodesic flow on G. The measure p, on M which is invariant under 
the geodesic flow Ut is the Haar measure dg of the group G. This measure 
is invariant under right and left multiplication by elements of the group. 
The Ut invariant scalar product for functions on G is therefore given by 



(/1I/2) = / fx(g)f2(g)dg 
Jgec 

2 r 00 />2tt />2tt 

4tt 2 J J q J 



f2(<f>,T,ip)f2(<f>,T,ip) sinhTdTd(/)dip. 



Finally we notice that right multiplication by elements in the rotation 
subgroup H = {e K ^} only affects the direction of momentum and does 
not change the position of the particle. Therefore we can identify the 
hyperbolic plane N with the set of right cosets G/H. 
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We now turn to the free motion on general surfaces of constant nega- 
tive curvature which are constructed by taking a tessellation of the hyper- 
bolic plane G/H and identifying the tessellating shapes (the fundamental 
domains). Let T be the group of transformations mapping the tessellat- 
ing shapes to each other. For every 7 g T the points jgH and gH of 
the hyperbolic plane are identified. The points of the quotient surface 
formed under this identification are labeled by double cosets in T\G/H. 
When the directions of momentum at each of the points are included the 
constant energy surface M = T\G is formed. Since the free motion is 
given by the right shift by h t it is not affected by quoticnting on the left 
by the subgroup T, so on the constant energy surface T\G the geodesic 
flow is given by 

U t (Tg) = Tgh t (19) 



4 Decomposition of the evolution operator 
for the geodesic flows 

In this section we use the representation theory of G — SO(2, 1) to find 
a decomposition of the form <fjh for the resolvent for the geodesic flow on 
T\G. 

The right regular representation Tr(/i) of G on L 2 (T\G) is defined by 

T R {h)ip{Tg) = ip{Tgh) where <p e L 2 (T\G). 

It can be decomposed as a direct sum of irreducible unitary represen- 
tations T v which leads to a splitting of L2(V\G) into a direct sum of the 
spaces W(T V ) on which T y acts. 

L 2 (T\G) = © M(T y ) (20) 



As the geodesic flow U t defined in (|19|) amounts to right multiplication 
by the group element h t the evolution operator Ut coincides with Tu(ht) 
and leaves M(T V ) invariant leading to the decomposition 

U t =Y,T y {h t ). (21) 

y£Y 

Substituting ( pT| ) into (||) we obtain the resolvent for the geodesic 
flow: 

/>oo 

R-(z) = - e- zt U t dt=Y J R-{ z ) ( 22 ) 
where R y _(z) = - J °° e - zt Ty(h t ) dt. 
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The rest of this section is organized as follows. In section 4.1 



we 



discuss the unitary irreducible representations T v . In section 4.2 we study 



the decompositions ( |20| ) and relate it to the spectra of the Laplacian on 
the quotient surface T\G/H. In section |4.3| we calculate the integrals 



R v _{z) for each of these irreducible representations and in section 4.4 we 
combine all the results to find R-(z). 



4.1 Irreducible representations of 50(2,1) 

Let T be an arbitrary unitary representation of G = 50(2, 1) on the 
Hilbert space H(T) with scalar product ( ). The Casimir operator is 
defined on H(T) by 

fi(T) = L%(T) + L 2 Q (T) - L\{T) 

where Lx(T) is the Lie derivative of T in the direction X G g and 
P,Q,K G g arc as defined in ([l?]). This operator commutes with each of 
the T(g) [|l5| and therefore must be a scalar multiple of the identity on 
each of the irreducible representations T y . It is convenient to denote by 
T p the unitary irreducible representation on which this scalar is — \ — p 2 , 
i.e. for all v G H(T P ) we have: 

fi(T> = (~-p> (23) 

The following values of p are allowed |T^| : 

• Im p = and Re p > (the principal series) 

• Re p = and Imp G (0, h) (the complimentary series) 

• Re p = and Im peN each corresponding to a pair of inequivalent 
representations (the discrete series). 

In addition we have the 1-dimensional identity representation I for which 
fi(i> = for v G H(J). 

Under the action of the compact abelian subgroup H = {e Kt } the 
representation T p splits into one dimensional irreducible representations. 
Let \n) G H(T P ) be a vector in one such one-dimensional representation 
with 

T p {e Kt )\n) =e lnt \n) (24) 

It is a fact (lj| that in a representation T p there is at most one vector 
|n) for each value n and they can be normalized to form an orthonormal 
basis for the space H(T P ). We also use the notation \p,n) when the 
irreducible representation to which \n) belongs needs to be specified. For 
the representations of the principal and complimentary series the value of 
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n ranges over all the integers. For the discrete series, one of the pair has a 
basis consisting of \n) where n > Imp and the other where n < Imp Q. 
A vector \<p) in H{T P ) may be expanded in the basis {\n}} as 

|¥>>=$>|P>|n> (25) 

n 

In the basis {|^)} the matrix elements for the representation T p are given 

by El 

(m\TP(g)\n) = e"^e^B l mn (co S hr) 
where I = ^ + ip and (cosh r) are the Jacobi functions. 

4.2 Decomposition of L 2 (T\G) 

We now consider the decomposition of L/2(T\G) into irreducible repre- 
sentations T p for a particular quotient surface Y\G/H. 

Suppose that the decomposition ( p0| ) of L2(r\G) has N p copies of the 
irreducible representation T p , labeled by s € {1, A^}: 

L 2 (T\G) = ® P (®f*MT p ' s ))- (26) 



Each H(T p;s ) has a basis {|p, n;s)} as described in section (4.1). From 
(p3) and (p3) we deduce that the space spanned by the vectors {|p, n; s) where s 



1, Np} is the intersection of the eigenspaces of Q(Tr) and Lk(Tr) with 
the eigenvalues — | — p 2 and m respectively. 

n(T R )|p,n;s) = (-i - p 2 )|p,n;s) 

ife(Tfl)|p,n.;s) = in|p,n;s) (27) 

The functions |p, 0;s) are invariant under multiplication on the right 
by elements of H and therefore can be viewed as functions in L 2 (T\G / H) . 
On this space the Casimir operator reduces to the Laplacian A hence 
the values of p occuring in ( p6| ) correspond to the eigenvalues e(p) — 
— j — p 2 of the Laplacian A on the quotient surface T\G/H. N p is the 
multiplicity of the eigenvalue e(p) and \p, 0; s) are its eigenf unctions. It 
is known that the spectrum Ar of A on a compact surface T\G/H is 
discrete and is bounded above by 0. It also contains the eigenvalue 
which corresponds to the the constant eigenfunction. Hence ( p6| ) takes 
the form: 

L 2 (r\G)= ® (®^iH(T" ;s )). (28) 
-i- p 2 eAr 

A differential equation for the eigenfunctions |p, n; s) is obtained by 
unwrapping a function in _L2(r\G) to give a periodic function on the 
whole group x '■ G — > C which obeys 

x(jg) = x(g) for all 7 e r. 
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Using Euler's coordinates ( |18[ ) on G in (]27|) we deduce that the unwrap- 
ping of the function \p,n;s) has the form Xs'"( T > 9)e m ^ where Xs' n ( T ^) 
obeys the second order differential equation: 

1 d sinhr^ + * (* - WnooBhr^ - n 2 ")) tf" = (-7 - p 2 )*^ 



sinhrdr dr sinh 2 t V<9 2 </> ^ // 4 

(29) 

4.3 Spectral decompositions for the irreducible rep- 
resentations 

We now find the correlation functions of the resolvents R p _{z) for each of 
the irreducible representations T p and obtain a spectral decomposition 
for the operator T p (h t ). 
Expanding 



Jo 

in the basis {\n)} we get 

/>OC 

(t\B£(z)\<p) = I e-^^\m)(m\T"(h t )\n)(n\^)dt (30) 

The matrix element (m\T p (h t )\n) can be written as a sum over the ex- 
ponentials e xt where 

A € {/,/- 1,/- 2...-JJ- 1,7-2,...} = A p (31) 

and 

(m\T"(h t )\n) = B l mn (cosht) = £ Cj A „„ e At (32) 
Substituting (|3|) into the resolvent (pp[) we obtain 



/>00 

<£|J£(*)Ip> = X! / e (A - z)t £<£M<i>l^ 



AeA p 

(CI TO ) C mn(«l^> 



EE' 



A6A p m,ra 

= £^>. (33) 

AGA P 

The operators .fiT^ determining the residue at the pole A £ A p are given 

by 



K x = J2\™)cL n (n\. 
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The matrix elements c„ n may be split into a product (appendix |Bj) c„ ln = 
a mPn which enables us to write (Q) in the form 

muzM = e - — 

AeAp 

= 2. x-z (34) 

AeA p 

where A runs through the set A p (^Tj) and 

|/A>=$>mH > I/-A>=E 6 ™H 
m m 

are eigenvectors of T p (ht): 

T"(h t )\f v )=e^\f v ). 

The coefficients and 6^ where X — 1,1 are found explicitly in appendix 
|b| along with a procedure for generating the eigenfunctionals for the other 
values of A e A p . 

The eigenfunctionals \f\) and |/_^) are linear functionals on the 
space § p C H(T P ) of test vectors \ip) = £ n c„|/9,n) e H(T P ) where 
the coefficients c n — > as |n| — > oo faster than any power of n i.e. 
lim^i^oo c n n k = for all k e N. By comparing (|I]) with (Q) we arrive 
at the spectral decomposition for T p (ht) 

T"{h t )= J2 e X Vx){f-x\- ( 35 ) 
AeAp 

This decomposition converges absolutely for t > for test vectors in the 
dense subspace T p c § p of vector; 
H(T P ) for some K. (appendix |cj). 



dense subspace T p C S p of vectors of the form \(p) — J2n=~K c n\p> n ) S 



4.4 Spectral decompositions of the evolution opera- 
tors for the geodesic flows 



On combining (35) with ( pq ) we obtain the central result of this paper: 
a spectral decomposition for C7 t : 

N p 

U t =T R {h t )= E E E eAt |/A ;S X/!xl (36) 

-i-p 2 eAr AeAp s=l 



where Ar is the spectrum of the Laplacian on T\G/H containing eigen- 
values — \~p 2 with multiplicities N p and where A p given by equation 
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( |3l| ) and also shown in figure [l]. The functional 

n 
n 

are the eigenfunctionals of U t : 

The functions \p, n; s), s = 1, ...N p form an orthogonal basis in the space 
of solutions of the linear differential equation (p9[). The coefficients are 
obtainable by expanding B l mn (cosh t) in powers of e xt (see equation |32| ) 
and are given explicitly for A = 1,1 in (p5|). 

The eigenfunctionals \f^' s ) belong to the space C°°(r\G)* and act 
on test functions in © § p = C°°(T\G). The decomposition (pi) 

-^-p 2 eAr 

converges for functions in a dense subspace 

N„ K 

T = {\<f) = £ £ £ c n , p . iS \ P> n;s) for some K}CC°°(T\G) 

(37) 

of test functions whose unwrappings have only a finite number of Fourier 
components in the angle ip (appendix |^) . 

The eigenvalues entering ([36]) could also be obtained by comparing 
classical orbit expansions for the traces of the resolvents of L and A |17j . 



5 Properties of the evolution operators for 
the geodesic flows 

5.1 Decay of correlations 

At long times only the leading terms with k = in each of the A p in the 
spectral decomposition for U t d3rj| ) remain significant (see section 2.4). 
Therefore as t oo, for f, ip e T c C°°(r\G) 

-I-p2 S A r s 

Separating out the eigenvalue with A = for a compact quotient surface 
where the dynamics is mixing we find that 
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The approach to equilibrium is governed by the eigenvalue(s) e^~5 ±J p)* 
in this expansion with the largest modulus (see figure [l]) . There are two 
possibile cases distinguished by the smallest non-zero eigenvalue — eo of 
-A. 

1. If — eo < j then the slowest decaying term gives the decay rate 
e A °* where Ao = — | + »J\ + eo. For the non-constant part of the 
correlation function (£\Ut\p) between £, p <E C oa (T\G) we have 

m t \p)-m(l\p)=0(e^) (39) 

2. If —cq > \ then all the values of p in (|38|) are real and we get 



A = -i 



4 

2 



(e|W-(e|l)(l|^ = 0(e^) (40) 



By applying the result ((TJ) to the geodesic flows we find that for 
arbitrary £,ip G C°°(M) 

mu t \v)-(a\i){i\p)\<ce^ 

This is a classical result proving the exponential decay of correlations for 
geodesic flows on constant negative curvature surfaces (see e.g. jl5) for a 
detailed discussion). Also from (|lj) we get the other terms in the long 
time asymptote: 

m t \p) - <eii)<%) 

N p N 

< C(N)e^ ^ N+1 ^. (41) 



5.2 Relation to the Laplacian 

In addition to refining the long time asymptotic form of the correlation 
function ( [i"l| ) the full expansion ( |36| ) provides a more detailed description 
of the evolution of distribution functions. Consider a smooth function 
^ on r\G which is independent of the momentum coordinate %j) (hence 
belongs to the space T) . As it evolves under U t this function will acquire 
some angular dependence but we will only be interested in the density 
on the configuration space i.e. the projection ^ t — PUt^ where P is the 
projector T\G T\G/H: 

1 f 2 * 
p X = 7T- / X dip 
27r Jo 
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Figure 1: The poles A — — i — ip — k (crosses) of the resolvent give 
the eigenvalues e xt of the evolution operator appearing in decomposition 
(|3(f). These poles are related to the eigenvalues — e = \+ p 2 of— A (open 
circles). The rate of decay e A °* at long times is determined by the pole 
Aq with the least negative real part. 



Expanding f as a linear combination of the eigenfunctions \p, 0; s) of 
the Laplacian for the quotient surface 

1*)= E £<P,<M*>|p,0;*> 

-i-p 2 £A r s 

and using 

U t \p,0;s) = TP(h t )\p,0;s) 

= ^2\p,n;s)(p,n;s\T p (h t )\p,0;s) 

n 
n 

combined with the following expressions for the matrix elements of the 
projection operator 

P\p,n-s) = n^O 
P\p,n;s) = \p,0',s) 

we obtain for 

\* t ) = Ptj t \*) = ]T Soo l+V (coshi)^(p,0, S |*)|p,0; S ) (42) 

Using the symbolic notation 

^-A- ||p,0;s) =p\p,Q;s) 
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we rewrite (|42|) as 



PU t \*) = bJ +1 ^ A ^(coshi)|*) (43) 



which relates the Laplacian A on the surfaces of constant negative cur- 
vature to the classical dynamics in the configuration space. 



6 Conclusions 

6.1 Summary of results 

Decompositions for the evolution operators Ut of a general mixing dy- 
namical system on a manifold M were constructed following Ruelle's pre- 
scription Q using the analytic continuation of the resolvent (^) : 

/>oo 

F^(z) = - e- zT m T \^)dT. 
Jo 

It was proved that when conditions ^ are met we have the follow- 
ing decompostions of Ut'. 

(£\U t \<p) = eM E^aI^I/a) convergent for t > 

A k 

{£\Ut\<p) = E e ~* E (/* I^Xel convergent for t < 

A fc 

where |/*) are eigenfunctionals of f/t 

Ut\f%) = e*\f*) 

belonging to the space C°°(M)*. They can be obtained from the residues 
of F^ tV> {z) which are positioned at A and have the form 



2™Res(A,F 5 , v ) = ^|/ A fc )(/^>). 



On applying this method to the evolution operator Ut for the free 
motion on a compact surface of constant negative curvature T\G/H we 
obtained the decomposition ( |36| ) which is convergent for £ and ip in a dense 
subspace T C C°°(T\G) of functions whose unwrappings onto functions 
of G have only a finite number of Fourier harmonics in the Euler angle tp 



(see equations 18, |37|). From this decomposition we obtained a refinement 
for the rate of decay of correlations Q) for these systems. 

We also found that the projection of the evolution operator PUt from 
T\G to T\G/H is related to the Laplacian A on the surface T\G/H by: 

Put = s -HV-A-i (cosht) 
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6.2 Consequences for the ballistic a-model 

Finally we discuss some consequences of the spectral decomposition (||) 
for the ballistic a- model [Q, Without diving into (the still controver- 
sial) issue of the cr-model derivation and its region of validity we quote 
below the result for the effective action (see [Q. |[ [L8| for a detailed dis- 
cussion) : 

S = JdfitrAW- 1 {L + i(u; + iO)A}w, (44) 

where W(x) belongs to some (super)-group and A is a particular matrix 
in this group obeying A 2 = 1 . The detailed structure of the target space 
is somewhat involved and for the purpose of our discussion it suffices to 
consider a toy model with W £ SU(2) and A = diag(l, —1). Then the 
target space is a two-dimensional sphere S 2 = SU(2)/U{1) parametrized 
by the matrices Q — W~ 1 AW. The action ( |44| ) does not depend on the 
parametrization of S 2 in terms of W. Indeed, introducing the two-form 
Q*V on the energy shell M obtained by pulling back the invariant volume 
V on S 2 by the function Q(x) : M — > S 2 we get for the action: 



S= (Q*V) A pdq + i(u + iO) tr(AQ)d/x, 

J M 

where pdq is the antiderivative of the simplectic structure dpAdq. If Q» V 
is exact (and it is certainly closed) the first term does not depend on the 
choice of the antiderivative but only on the simplectic structure dp A dq. 

It is universally believed that some sort of regularization should sup- 
plement the effective action (|4|). It was conjectured Jj| that in the limit of 
a vanishing regulator the eigenvalues of the regularized Liouvillian opera- 
tor approach the Rucllc resonances (which were refered to as "eigenvalues 
of the Perron Frobenius operator" in S). 

We would like to point out, that this conjecture must be further clari- 
fied due to the existence of the two inequivalent sets of Ruelle resonances. 
One set of resonances is in the left half of the complex plane, while the 
other is in the right. These two sets originate from the different branches 
of the resolvent given by the two integral representations (equations ^| 
and [l5] respectively). As a result there exists two non-equivalent regu- 
larizations of the operator L. Denoting by L reg the regularization of L 
with eigenvalues close to the Ruelle resonances in the left half plane and 
observing that the operator —L+ eg has the eigenvalues close to the Ruelle 
resonances in the right half plane we suggest that the two regularizations 
of L are L reg and ~L^ eg . 

We suggest the following structure for the regularized operator in the 
target space: 
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which ensures the convergence of the action ([44]). Expanding the ac- 
tion (Q) near Q — A as W = 1 + ^ ^_ ™ j we verify that the choice 

( |45| ) ensures that the quadratic part of the action S 2 S is non-postively 
defined: 

^ -/*<•-) (A 

Note, that a similar structure for the regularized action appeared in the 
model with diffusive scattering p9[. 
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A Mixing dynamical systems 

A dynamical system (M, Ut, fi) is mixing if for any two sets A, B C M 
lim n(U t A nfl)= K A )K B ) (46) 

i — >oo 

For a mixing system the only square integrable eigenfunctions of U t 
are the constant functions which have the eigenvalue 1. The proof |2(| of 
this result is reproduced below. 

We write (|6|) in terms of the characteristic functions \A and \b of 
the sets A and B (The characteristic function xc of a set C takes on the 
value xc{ x ) = 1 if % £ C an d otherwise). 

lim (& X a|xb) = (xa|1)(1|xb) (47) 

t — >oo 

Since the space of linear combinations of characteristic functions is dense 
in the space of square integrable functions ( [f7|) must also be true for 

f,geL 2 (M) 

]im(U t f\g) = (f\l)(l\g) (48) 

t — >oo 

Let / be an eigenfunction of Ut with eigenvalue e xt and 5 = 1 then jig ) 
becomes 

e At (/|l) = </|l) 
therefore the eigenvalue must be 1. 
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B Coefficients for the eigenfunctionals 

We now study the matrix elements c^ ln = (m\K\\n) of the operator K\ 
appearing in the residues of the resolvent (33) 

,m rP , z ^ s = (£\K\\<p) = y- (g|m)c* n (re|y>) 
r~f A - z /-f ^ A - z 

AeA p AeAp m.n 

Consider the following operators on H(T P ) 

L = L P (TP) 
B_ = L K (TP)+L Q (TP) 
B + = L K {TP)-L Q {TP). 

Note that T p (h t ) = e Lt and we have the commutation relations 

[L,B+] = B+ 
[L,B-] = -B_ 
[B-,B+] = 2L (49) 

From the commutation relations we have that B + L = (L+l)B + there- 
fore B + e Lt = e^ L+1 ^B+ i.e. B + TP(h t ) = e t TP(h t )B + . By considering 
(£\B + R p _\(p) and using this result we obtain 

B+k x = K X+1 B+ (50) 

Hence B + sends the image of K\ into the image of K\+\. In particular 
B + Ki = = B + Kj = since there are no eigenvalues I + 1 and I + 1 in 
the set A p . 

Consider an eigenfunctional |/) in the image of K\ so that it satisfies 
TP(h t )\f) = e xt \f) and therefore 

£|/>=A|/>. 

Suppose also that it is annihilated by B+: 

B+\f)=0 

so that using (^) we get 

(n\B+\f) = 
(n\B + B_\f) = -2{n\L\f) = -2A(n|/) 

for all n E Z. 

Representing |/) by the linear combination 



I/) = 51 a n \n) where a n = (n\f) 
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and using the explicit expression for the matrix elements of B + and £?_ 


-2i(n\B + = -(l + n)(n-l\-2n(n\ + (l-n)(n+l\ 
-2i(n\B- = (I + n)(n - l\ - 2n(n\ - (I - n)(n + l\ (51) 

we get the system of linear equations 

— (I + n)a n -x — 2na n + (I — n)a n +i = 

-(l + n- l)(l + n)a„_ 2 - 2(1 + n)a n - X + 2(1(1 + 1) + n 2 )a n 

- 2(1 - n)a n+ i - (l-n - l)(l - rz.)a„_ 2 = -4Aa„ 

These equations have only the trivial solution a n = unless A — I — 
— ^ + ip or A =7 = — | — ip in which case they have the unique non-trivial 
solutions 

oL - 1 



(l + n)\(l-n)\ 



(~ir (52) 



respectively. Hence the images of Ki and Kj are one-dimensional. 

Since B + sends the image of K\ into the image of K\+\ and its 
restriction on the image of K\ has zero kernel for A ^ 1,1 we deduce that 
the images of Ki-k and Kj_ k for all k £ N are also one-dimensional. The 
operators K\ where A e A p may therefore be written as 

= |/XX/_I| = E M a m&£H = E N^H 



and their matrix elements c^ nn split into a product c^ nn — a^b* where 

\h) = E a -i™> 

m 

l/-A> = E 5 » ( 53 ) 

n 

as stated in the main text. 

Combining (52) with (53) we see that a l n = a l n and a l n = a~ n . This 
can be confirmed from equation d32| ) and the expansion of B l m0 (cosh t) at 
large t |l| : 

/2 l (V) 2 1 

BL,(cOBhi) = (—^-75 7771 T7 e " 

mtn ; I Vif (Z + m)!(J-m)! 



2\l\) 



(-l) m e") 



(l + 0(e-*)). (54) 
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Using (d4) and the symmetry relation p 



j . . a + n)!(/-n)! , , 

B l mn (cosht) = -B l nm coshi 

(t + rnj!(( — my. 



we can find the coefficients b l n and 6^: 

£ = ^L(/!) 2 G + n)!(/-n)! 

V 7r 

To summarize the functional |/a), |/_t) involved in the decomposi- 
tion of K\ for A = 1,1 have the following form 

I/') - E (j + _„).!"> 

i/t> = £(-i)» 

= £^=(Z!) 2 |n)and 

— V 71 " 
n v 

\f-l) = ^^=(n) 2 (/ + n)!G-n)!|n) (55) 

Given the eigenfunctionals |/z) and the eigenfunctionals |/x) for 
the other values of A e A p may be obtained by successive applications of 
the operator £?_ : 

\f x _ k )=B k \f x ) (56) 

due to the relation B^K\ = K\-\B_ which is proved analogously to 
©■ 

Similarly, the projectors can be obtained by applying the operator B + 
to |/_ T ) and |/_,): 

\f. J+k )=B k + \f_j). (57) 
From (O) and the relations (KqM) we find that 



a l - k = (n|/i_ fc ) = (n\B k _\fi) = 0{n k ) 

a\7 k = (n\f 1 _ k ) = (n\B k _\f T ) = 0(n k ) 

b l ~ k = (n\f_- l+k ) = (n\B k + \f_j) = 0(n k ) 

b J - k = (n\f_ l+k ) = (n\BX\f_ l )=0(n k ) (58) 
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C Convergence of the spectral decomposi- 
tions 

We now present the necessary condition for the convergence of the spec- 
tral decomposition (|35l): 

W(h t )\<p) = e xt {f_ T \ip){i\h). (59) 
AeA p 

In particular we show it is convergent when £ and tp belong to the dense 
subspace T p of S p C H(T P ) where 

K 

Y p = {\tp)= c n\p,n) for some K} (60) 

n=-K 

First we show that each term e xt (f_j\ip}(^\f\) in fl59| ) is defined for £ 
and </? in the subspace § p : 

§p = {\<p) = V c„ |/9, n) where lim c„n« = for all q e N}. (61) 

— ' In — >oc 

n 

Note that § p is invariant under T p (h t ). Let £ and </? be given by 

\<f) = ^2c n \p,n) 
n 

\0 = Y,d n \p,n) (62) 

n 

Using (f_j\m){n\fx) = a^b* we get 

|e Xt (/_ 5 |v)<e|A>| < E l^iKcM (63) 

The eigenvalue A is of the form I — k or 7— k so by ( |5E} ) = 0{m k ) and 
6* = 0(n k ). Hence 

|e At (/_ x |^)<£|/A)| < Cle^l^l^l^lW 

tci n 

which converges due to (|6l]). Therefore each term is defined for £, <y9 6 §> p 
and the eigenfunctionals |/„) belong to the space § p *. 

Now we show that the sum ( |59| ) converges for £,tp S T p . Since the 
sums (|6^) for £, ^ have only a finite number of terms it suffices to prove 
convergence for \ip) — \n) and |£) = \m). In this case fl5S| ) reduces to (^) 
which converges absolutely for t > 0. 
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Now we turn to the convergence of the spectral decomposition of the 
evolution operator U t ( |36| ) 

(t\u*\<p)= E E E ^W)(f%\<P) ( 64 ) 

-±-p 2 GA r 8=1 AeA p 

By the above considerations the terms in this series e At (£ |/X' s )(/ P y|y) 
are defined for f,<p € ffi § p = C°°(r\G). 

We will show that expansion ( p^ ) is convergent when £, </? belong to 
the subspace T c C°°(T\G): 

N p K 

T = {\f) = E EE c n,p;s\p,n;s) for some K} (65) 

-i-p 2 GA r s =1 n=-K 

It is sufficient to determine the convergence for 

w) = E E c ^ s i^ n;s ^ 

-i-p 2 GAr s=1 

10 = E E d p; s ^' m;s ) 

-i-p 2 GAr 5=1 

for which decomposition ( p4[ ) is 

iv P 

(Cl&M = E E E e At c p , fl a*d p , fl 6* (66) 
-i-p 2 SAr 8=1 A e A P 

Let tp(T,9,i/j) = (p n (T,8)e m ^ be the unwrapping of \ip). The function 
(,o„ on T\G/H has the norm 

||Vn|| 2 = — / \Pn( T , S)\ 2 smh rdrdO , F is a fundamental domain 

27T Jp 

We may expand ip n in the eigenfunctions Xn S (the unwrappings of 
1/5, n; s)) of the operator (29): 

<p n (r,o)= E E c ^« s - ( 6? ) 

-i-p 2 GA r »=1 



hence 



irf = E Ei^i 2 - 

-i-p 2 SA r 8=1 
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Similarly 

N p 

\\u\ 2 = E Ei^i 2 - 

-i-p 2 eA r s=1 

Using the fact that 

N P 

E Ei c p.' d ".«i<ii^iiii^ii 



,e At l 



and that ( |32| ) converges absolutely we see from (|66|) that 



mu t \<p)\ = E Ei c p^i E ^ 

-i-p 2 6Ar s=1 A eA p 

< c E Ei c p^i <cy ii^iiii^ii 

-i-p 2 GA r «=1 

is bounded. 

Hence taking £, 99 € T ensures the convergence of the spectral decom- 



position (64). Since the unwrapping of \p, n;s) has the form Xn" the 
unwrappings of functions in T will have only a finite number of Fourier 
components in ip. 
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